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1.  Introduction 

We  address  the  problem  of  computationally  efficient  representation  of  polyhedra. 
Many  algorithms  for  manipulating  polyhedra  are  naturally  decomposed  in  two  steps: 
the  localization  of  the  faces  that  intersect  a  given  field  of  interest,  and  the  subsequent 
processing  of  these  faces.  Here  are  some  practical  examples: 

1.  Clipping:  Localize  all  the  faces  of  a  polyhedron  that  intersect  a  viewing  pyramid, 
then  process  them  to  determine  which  ones  are  visible. 

2.  Set  Operations  between  solids:  Localize  all  the  intersecting  faces  of  two  polyhedra, 
then  jtrocess  them  to  find  the  intersection  polygons  of  the  two  polyheclra,  and  to 
determine  whether  or  not  all  the  other  faces  are  part  of  the  resulting  solid. 

3.  Ray  Casting:  Localize  all  the  faces  of  a  polyhedron  intersected  by  a  semi-infinite 
straight  line  (the  light  ray),  and  process  these  faces  to  find  the  first  intersection,  and 
compute  the  display  parameters. 

In  the  case  of  polyhedra  with  many  faces  (typically  more  than  one  thousand),  the 
localization  step  is  the  most  expensive.  This  is  more  generally  true  of  algorithms  that 
manipulate  a  large  number  of  objects.  The  basic  method  proposed  in  the  litterature 
to  speed  up  this  step  is  to  include  each  object  in  a  simple  enclosing  box,  thei  build  a 
hierarchy  of  these  boxes,  and  use  it  to  direct  the  search.  This  method  has  beer  used  for 
the  Clipping  and  Ray  Casting  problems  by  Clark  [1976],  Rubin  and  Whitted  [1980],  and 
more  recently  by  Weghorst,  Hooper,  and  Greenberg  [1985].  It  has  also  been  applied  to 
algorithms  for  Set  Operations,  essentially  in  the  Constructive  Solid  Geometry  approach 
(Requicha  and  Voelcker  [1982]). 

All  these  algorithms  are  "object”  based:  low  level  boxes  in  the  hierarchy  usually 
represent  either  complete  polyhedra,  or  primitive  solids.  They  are  not  well  suited  to 
the  manipulation  of  a  set  of  polyhedra  each  having  thousands  of  faces.  Such  polyhedra 
are  of  practical  importance,  in  particular  they  can  be  models  of  real,  digitized  objects 
that  no  simple  CAD  system  can  handle.  Examples  are  organs  from  CT  scanners  (Artsy, 
Frieder  and  Hermann  [1978]),  or  objects  measured  using  a  laser  rangefinder  (Faugeras 
and  Pauchon  [1983]). 

Two  recent  methods,  one  for  Set  Operations  between  solids  (Mantyla  and  Tamminen 
[1983]),  and  the  other  for  Ray  Casting  of  fractal  surfaces  (Kajiya  [1983])  tave  been 
designed  to  deal  explicitly  with  such  polyhedra.  They  use  “face”  hierarchies  that  enclose 
the  faces  themselves  in  boxes. 

We  present  a  new  method,  called  Prism  Tree ,  to  represent  polyhedra  with  many 
faces.  This  is  also  a  “face”  representation,  that  describes  polyhedra  as  a  ternary  tree 
of  prisms  as  enclosing  boxes  (Figure  1).  It  generalises  the  Strip  TVce  represer  tation  for 
planar  curves  (Ballard  [1981])  and  derives  from  a  polyhedral  approximation  algorithm 
(Faugcrai  et  al.  [1984])  that  is  described  in  Section  2.  However,  Prism  Trees  lave  much 
wider  applications,  in  particular  they  are  information  preserving,  which  is  no:,  the  case 
of  the  approximation  algorithm. 

The  representation  is  intrinsic  to  the  surface,  as  opposed  to  the  EXCELL  structure 
used  by  Mantyla  and  Tamminen  [1983],  or  to  the  Octree  representation  for  volumes 


Figure  1.  The  Prim  Tree  representation  for  paqrhedra.  At  each  level,  the  surface  ST  associated  with  a 
node  is  represented  by  an  enclosing  prism  (dotted  lines)  The  surfaces  associated  with  the  sons  of  a  node 
partition  the  surface  associated  to  their  father.  The  link  between  a  node  and  the  associated  ST  is  only  made 
explicit  at  the  leaf  level  (solid  fines).  As  these  surfaces  partition  the  original  polyhedron,  the  representation 
is  information  presening. _ _ 

(Jackins  and  Tanimoto  [1980],  Meagher  [1982],  Iftikhar  [1981]),  that  are  attached  to  a 
particular  reference  frame.  It  is  therefore  invariant  under  rigid  transformations.  We  show 
that  Prism  Drees  can  be  used  to  speed  up  the  localization  step  of  Set  Operations  and  Ray 
Casting  algorithms  for  polyhedra  with  many  faces.  They  unify  the  approaches  of  the 
Mantyla  and  Tamminen  (although  the  data  structures  are  very  different)  and  Kajiya 
algorithms,  and  in  general  demonstrate  an  improvement  over  these  methods. 

We  first  sketch  the  underlying  polyhedral  approximation  algorithm  (Section  2).  Then 
we  define  the  Prism  Trtt  structure,  and  give  some  properties  (Section  3).  Section  4  de¬ 
scribes  a  set  operations  algorithms  for  Prism  TYees,  and  Section  5  generalises  to  the  Prism 
Trtt  structure  a  Ray  Casting  algorithm  proposed  recently  by  Kajiya  [1983]  for  fractal  sur¬ 
faces.  The  crucial  localization  step  has  been  fully  implemented.  The  implementation  of 
the  processing  step  is  under  way.  Examples,  featuring  set  operations  between  synthetic 
solids  and  Ray  Casting  display  of  both  synthetic  and  real  objects,  are  given.  Detailed 


procedural  versions  of  the  algorithms  can  be  fonnd  in  the  appendix. 

8.  Pd^hedral  Approximation 

la  this  section,  we  describe  the  polyhedral  approximation  algorithm  (Fauge.as  et  al. 
[1984]).  As  we  will  see  in  the  sequel,  the  Prism  Tree  employs  a  variant  of  this  algorithm. 
The  diaia  consists  of  a  polyhedron  of  genus  0  (ie  without  holes,  we  will  see  the  reason 
of  this  restriction  below),  whose  faces  are  not  necessarily  triangular.  The  polyhedron  is 
described  by  an  object  graph  OG  whose  nodes  are  the  vertices  of  the  polyhedron,  and 
arcs  are  the  polyhedral  edges  joining  them.  The  algorithm  generalises  to  3D  space  a 
recursive  polygonal  approximation  algorithm  (Duda  and  Hart  [1973]).  Using  a  breadth 
first  approach,  it  approximates  the  initial  polyhedron  by  a  polyhedron  with  triangular 
faces  Tit. 


Figure  2.  The  split  step.  The  triangle  T  =  PQR  is  replaced  by  the  three  triangles  T\  =  PQM.Tt  = 
QRM,Ti  -  RPM  that  are  closer  to  the  original  surface. _ _ 


At  each  level  of  the  recursion,  we  associate  to  each  triangle  T  a  set  of  attributes:  a  part 
ST  of  the  surface  to  approximate,  the  three  vertices  P,Q,R  of  T,  the  error  measured 
by  the  maximum  distance  between  the  points  of  ST  and  the  triangle  plane,  anti  ;he  point 
M  where  this  maximum  is  reached.  The  triangulation  is  represented  by  an  adjacency 
graph  AG:  the  nodes  are  the  triangles  themselves,  and  two  nodes  are  linked  by  an  arc  iff 
they  share  a  common  edge. 

Starting  with  a  rough  approximation  (described  below)  of  the  surface  ,  the  algorithm 
loops  over  the  following  steps  until  the  error  associated  to  each  triangle  is  less  than  a 
given  threshold  e. 


Prim  Tree* 


Split  Stqn  For  each  triangle  T  such  that  ey  >  e,  do:  (Figure  2) 

I.  Let  Tx,Ti,T%  be  the  three  triangles  PQM ,  QRM  and  RPM,  and  STi,STi,STz  the 
mirriitnrt  surfaces  (obtained  with  a  method  to  be  described),  compute  the  corre¬ 
sponding  errors  ey,,  ey,,  er,,  and  the  points  where  they  are  reached 


2.  Replace  T  by  the  X<s  in  the  adjacency  graph  (in  particular,  for  each  previously  neigh¬ 
bor  of  T,  replace  T  in  its  list  of  neighbors  by  the  coresponding  Tj). 

If  we  only  use  this  step,  old  edges  are  never  removed,  even  if  they  are  a  very  bad  approx¬ 
imation  of  the  surface.  This  is  the  reason  why  we  use  the  following  step. 


Figure  3.  The  Adjustment  step.  The  triangle*  T*  and  T*  are  replaced  respectively  by  T*  and  Tf,  and  T, 


Adjustment  Step:  For  each  pair  of  neighboring  triangles  Ta  =  PaQR  and  T*  =  P*QR 
that  have  been  created  at  the  previous  step  and  whose  associated  error  is  greater  than  e, 
eliminate  the  edge  QR  (Figure  3): 

1.  Find  in  ST *  U  ST*  the  point  M  that  is  the  closest  to  the  bissector  plane  of  T*  and 
1 *,  and  lies  at  the  maximum  distance  from  QR. 

2.  Let  7^  and  Tf  be  the  triangles  P°QM  and  PaRM}  compute  as  before  the  associated 
surfaces  STf,  STf,  the  errors  and  f j,  and  the  points  M°  and  M%.  Compute  the 
same  way  the  triangles  T*  and  T *  and  their  attributes. 

3.  Replace  T*  and  T*  by  the  four  new  triangles  in  AG. 
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PolyWdrnl  Approximation 

We  have  outlined  the  algorithm,  we  now  initialise  the  adjacency  graph,  and  d  ifine  the 
STn  recursively  . 

Initialisation  of  AG:  Choose  three  points  P,  Q,  R  on  the  surface  (Figure  4).  We  can 
compote  the  shortest  cycle  of  OG  that  contains  these  points  and  is  the  closes!  to  their 
plane.  As  the  surface  is  of  genus  0,  this  cycle  cuts  it  into  two  connected  compon  ents  STi 
and  STi  (Giblin  [1977]).  The  initial  graph  is  composed  of  two  nodes  correspcnding  to 
the  same  triangle  PQR,  but  whose  associated  surfaces  are  STi  and  STj.  The  t  vo  nodes 
are  linked  by  three  arcs  corresponding  to  the  three  edges  of  the  triangle.  Notice  though 
that  any  other  initial  triangulation  could  have  been  used.  In  particular,  trian  julations 
based  on  .he  smooth  patches  between  surface  discontinuities  (as  in  the  Intrinsic  Patches 
representations  of  Brady,  Ponce,  Yuille,  and  Asada  [1984])  could  be  used  to  yield  better 
approxim  itions  of  the  original  surface. 


(PQ1 


Triangle:  T 
Error  E  t 
Split  Point:  A/j 
Surface:  ST  j 

f  m 

Triangle:  T 
Error:  E  2 
Split  Point:M  2 
Surface:  ST  2 


Figure  4.  The  initial  adjacency  graph.  The  three  paths  [PQ],  [QR]  and  [RP]  define  a  cycle  tut  cuts  the 
surface  into  two  part?  STi  and  STj.  They  define  the  initial  graph. _ _ _ 


Definition  of  the  STs:  We  consider  only  the  split  step  (the  case  of  the  adjustment 
step  is  quite  similar  and  will  be  omitted),  and  use  once  again  the  fact  that  a  cycle  of  a 
surface  of  genus  0  cuts  this  surface  into  two  components.  Suppose  (Figure  5)  th  v.  we  have 
associated  to  a  triangle  T  a  surface  ST,  delimited  by  a  cycle  composed  of  the  ta  ee  paths 
[PQ],  [Qf?],  and  [RP\.  Let  P1.2,  P2.3,  P3.1  be  the  bissector  planes  of  the  pairs  (Ti.Tj), 
(T2,T3),  and  (T3,T,).  Find  a  path  [PM]  between  P  and  M  in  ST  as  the  shortest  path  as 
close  as  possible  to  P3,j.  Define  similarly  the  path  \QM\.  The  three  p<vths  \P.V  j,  [MQ], 
and  [ QP ]  define  a  cycie  \PMQ\  in  ST  that  cuts  it  into  two  connected  compcnents  ST\ 
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Figures.  The  wfoce  ST  i»  iplit  into  three  parts  by  the  paths  [PM],  [QM],  *nd[HM] 


and  ST.  We  associate  ST\  to  T\.  We  can  now  find  in  ST  a  path  [/iM],  that  splits  ST 
into  two  components  5T2  and  ST3,  that  we  associate  to  T3  and  T3. 

The  polyhedral  approximation  algorithm  (Faugeras  et  al.  [1984])  was  initially  pro¬ 
posed  in  the  context  of  Computer  Vision,  with  the  purpose  of  obtaining  from  an  original 
polyhedron  (typically  with  5000  faces)  a  reasonable  approximation  with  only  typically 
500  faces.  Tile  main  point  was  to  compress  the  information  as  much  as  possible,  while 
preserving  u  much  as  possible  the  shape  of  the  object.  Unfortunately  this  is  not  always 
possible,  in  particular  for  very  complicated  objects  with  important  concavities,  and  in 
practice  the  approximation  can  become  relatively  poor  (Ponce  1983]).  We  now  derive 
from  this  algorithm  the  Prism  Tree  structure,  and  show  that  these  problems  disappear, 
mainly  because  this  representation  is  information  preserving. 


3.  Prism  Tress 


3.1.  Definition 

The  polyhedral  approximation  algorithm  provides  us  naturally  with  a  hierarchical 
representation  of  the  surface:  the  successive  subdivisions  of  the  surface  yield  a  ternary 
tree  structure,  where  each  node  represents  a  triangle  and  the  associated  surface,  and 
stores  two  kinds  of  information,  geometric  and  structural. 

The  information  carried  by  the  STs  is  too  complex  for  being  directly  used  at  each 
level  of  the  tree,  so  the  “geometric”  part  of  a  non-leaf  node  is  a  simple  box,  easier  to 
manipulate.  However,  we  store  in  each  leaf  of  the  tree  a  description  of  the  .associated 
connected  component  ST  of  OG.  In  fact  this  description  will  not  only  conta  in  the  edge 
information  corresponding  to  OG,  but  also  a  description  of  the  faces  of  th2  original 
polyhedron  adjacent  to  ST.  This  is  important  as  it  is  in  fact  this  face  information  that 
will  be  used  in  the  processing  step  of  our  algorithms. 

Since  by  definition,  the  STs  partition  the  original  graph,  this  makes  the  rtp'vsentation 
information  preserving .  the  original  polyhedron  can  be  recovered  from  its  P'ism  Tree 
representation.  Consider  a  triangle  T,  and  the  bissector  planes  Pi, P‘2,  P3  of  T  and  its 
neighbors.  We  define  the  box  associated  with  ST  to  be  the  smallest  truncate  c  pyramid 
(called  prism)  with  three  faces  PP{  parallel  to  the  P{6,  and  the  two  remaining  ones,  TB\ 
and  TB2,  parallel  to  T,  that  encloses  ST  (Figure  6).f 

The  structural  information  of  each  node  consists  of  three  pointers,  noted  Son[t]> 
t  e  [1, 2, 3].  The  root  of  the  tree  points  to  the  two  half-surfaces  defined  by  the  initialisation 
step.  The  Sons  of  a  node  associated  to  a  split  triangle  point  to  its  three  sub- triangles, 
and  the  Sons  of  a  node  associated  to  an  adjusted  node  point  to  its  two  sub-triangles  (the 
third  pointer  is  nil).  We  will  refer  to  a  node  by  a  pointer  Pt  to  it.  Figure  7  shows  the 
Prism  Tree  model  of  a  sphere.  Notice  that  the  prisms  get  thinner  and  thinner  as  the 
resolution  increases. 

3.2.  Some  Properties 

We  now  give  fundamental  properties  which  relate  the  intersection  of  two  prisms  to 
the  intersection  of  the  underlying  surfaces.  The  first  one  is  obvious:  if  two  pi  isms  don’t 
intersect,  then  the  underlying  surfaces  don’t  intersect  either.  In  this  case,  wt  say  that 
the  associated  nodes  have  a  null  intersection. 

The  converse  proposition  is  evidently  false  in  general:  two  surfaces  can  hav;  an  empty 
intersection,  although  the  associated  prisms  intersect.  In  fact,  it  is  impossible  to  decide 
that  the  surfaces  associated  to  two  nodes  intersect  before  the  leaf  level,  as  the  STs 
themselves  are  only  represented  at  this  level.  So  we  say  that  two  nodes  that  are  not 

f  This  is  a  natural  way  to  define  the  boot  as,  by  definition,  ST  is  the  part  of  the  surface  th.it  lies  within 
the  Pi s  (or  more  precisely,  it  is  the  part  of  the  surface  delimited  by  a  cycle  which  is  the  intersection  of 
the  surface  and  the  three  planes),  even  though  in  the  general  case,  points  of  ST  ran  lie  outside  of  the  P,s 
(remark  however  that  if  ST  is  convex,  then  it  lies  entirely  inside  the  P,s).  The  two  remaining  faces  are 
used  to  “close”  the  box.  and  their  distance  corresponds  to  it- 


Prism  Trees 


Figure  6.  Definition  of  the  prism  associated  to  a  triangle.  The  upper  part  of  the  figure  shows  the  geometry 
of  the  biaaector  planes  PP,s  and  the  parallel  planes  TB{ s.  The  lower  part  shows  T  and  its  neighbors  T,s, 
and  the  constructed  enclosing  prism.  _ 


both  leaves  have  a  possible  intersection.  At  the  leaf  node,  we  can  decide  if  the  surfaces 
associated  to  two  leaves  Pt\  and  P<2  intersect  by  testing  directly  each  face  of  ST\  against 
each  face  of  572.  If  any  of  these  couples  intersect,  then  STi  and  ST2  intersect,  and  we 
say  that  Pti  and  Ptj  have  a  clear  intersection.  Otherwise,  we  say  again  that  they  have 
a  null  intersection. 

Let  ns  state  two  more  properties  that  are  obvious,  but  are  the  basis  of  all  the  algo- 


m 


rithms  we  have  developed  for  manipulating  Prism  TYees.  If  the  surface  associated  with 
a  node  intersects  an  object  (a  line,  another  surface..),  then  the  same  property  holds  for 
all  the  ancestors  of  the  node.  Conversely,  if  the  surface  associated  to  a  node  does  not 
intemect  an  object,  then  the  same  property  holds  for  all  the  descendants  of  this  node. 
This  will  allow  the  efficient  localization  of  intersections  by  pruning  the  tree  of  possibilities 
as  soon  as  possible. 


3.3.  A  lemma 

We  have  noted  in  the  previous  section  that  the  geometric  information  carried  by  a 
non-leaf  node  of  a  Priam  Tree  is  too  poor  to  directly  test  the  intersection  of  the  under¬ 
lying  surface  with  another  object.  We  now  give  a  lemma  that,  given  certain  additional 
conditions  of  regularity,  relates  the  intersection  of  a  straight  line  and  a  prism  to  the 
intersection  of  this  straight  line  and  the  surface  ST  associated  to  this  prism. 
Definition:  A  prism  is  said  to  be  regular  if  the  associated  planes  PPi,  PP2,  and  PP$ 
are  the  bissector  planes  Pi,  P2,  and  P3  themselves. 

The  regular  prism  notion  is  analogous  to  the  notion  of  regular  atrip  introduced  by  Ballard 
[1981]  for  Strip  Trees.  It  ensures  that  the  intersection  of  the  surface  ST  associated  to 
a  regular  prism  and  the  union  of  the  PP,-s  is  a  closed  curve.  It  is  then  easy  to  show 
the  following  lemma  (see  Faugeras  and  Ponce  [1985]  for  a  more  general  statement  and  a 
detailed  proof,  based  on  Jordan’s  theorem). 

Clear  Intersection  Lemma:  If  a  straight  line  intersects  both  the  planes  TB\  and  TBj 
of  a  regular  prism,  then  it  intersects  the  underlying  surface  ST  an  odd  number  of  times. 
As  it  implies  that  the  line  and  the  surface  intersect,  we  also  say  that  the  straight  line  and 
the  prism  have  a  clear  intersection  (Figure  8). 

This  lemma  could  be  used  for  deciding  if  a  point  is  inside  or  outside  a  polyhedron 
represented  by  a  Prism  Tree,  by  counting  the  number  of  clear  intersections  of  a  semi¬ 
infinite  straight  line  traced  from  this  point  with  the  nodes  of  the  tree,  using  a  method 
similar  to  the  one  developed  by  Kalay  [1982]  f. 

A  straight  line  that  has  a  clear  intersection  with  a  prism  intersects  at  least  one  time 
the  associated  ST.  We  have  implemented  this  property  in  our  Ray  Casting  algorithm 
(see  Section  5)  to  prone  as  soon  as  possible  the  list  of  the  triangles  that  may  intersect  a 
given  ray. 

3.4.  Soma  remarks 


Contrarily  to  the  polyhedral  approximation  algorithm,  the  Priam  Tree  is  proposed 
in  the  context  of  Geometric  Modelling,  and  it  is  used  only  for  localization  purpose.  In 
our  algorithms,  the  ultimate  processing  step  is  the  same  one  as  for  a  classical  polyhedral 
representation.  This  slightly  different  approach  has  some  important  consequences. 

First,  we  use  in  fact  a  variant  of  the  approximation  algorithm.  We  no  longer  stop 
dividing  a  node  when  the  error  gets  small,  but  when  the  site  (number  of  points)  of  the 
associated  ST  is  less  than  a  given  value  (6  points  per  leaf  in  the  examples  presented  in 


f  Mantyla  and  Tamminen  [1983],  using  their  Dcoc-EXCELL  data  structure,  have  actually  implemented 
a  similar  method  for  testing  the  inclusion  of  a  point  inside  a  polyhedron. 


m. 


T«  !■ 


Prism  Aces 


11 


Figure  8.  The  dear  intersection  lemma.  The  ray  intersects  ST  between  its  two  intersection*  with  the  7  i?,*. 


the  sequel):  so  the  depth  of  the  tree  representing  a  n  faces  polyhedron  is  0  (Log  (n)), 
and  the  number  of  faces  processed  is  proportional  to  the  number  of  leaves  found  during 
the  localization  step. 

Second,  as  the  representation  is  information  preserving,  and  as  the  Prism  Tree  is  used 
for  localization  only,  even  poor  approximations  of  an  object  give  exact  results.  The  only 
difference  is  that  the  tree  of  a  poorly  approximated  object  may  have  bigger  prisms  at  a 
given  level  in  the  tree  than  a  well  approximated  one,  and  that  the  tree  can  be  unbalanced. 

Note  also  that  the  Prism  Tree  is  intrinsic  to  the  surface  it  represents:  its  geometric 
features  are  not  constructed  according  to  a  particular  fra me,  but  to  surface  features.  This 
implies  that  the  structure  of  a  prism  tree  is  invariant  through  rigid  transformations. 

A  last,  but  important  remark:  the  approximation  algorithm  has  been  designed  for 
surfaces  of  genus  0,  which  is  an  unfortunate  restriction.  We  can  extend  our  algorithms 
and  analyses  to  more  general  surfaces  by  initially  dividing  them  in  components  of  genus 
0  (e.g.,  cut  a  torus  in  two  components  along  its  parabolic  lines),  and  considering  a  simple 
triangulation  of  these  components  as  the  initial  AG.  The  drawback  is  the  addition  of  a 
(possibly)  interactive  first  step  to  define  this  triangulation.  But  as  long  as  this  triangu- 
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laltea  is  simple  enough  (ie  the  number  of  initial  nodes  of  the  Priam  Tree  is  bounded  and 
«n»»ll  gygh),  it  does  not  change  the  overall  complexity  of  the  Priam  Tree  algorithms. 

The  price  can  be  worth  paying  for  including  very  complicated,  real  objects  in  a  high 
level  Graphics  environment.  This  is  quite  similar  to  the  approach  advocated  earlier,  that 
first  segments  the  surface  into  smooth  patches  (Brady,  Ponce,  Yuille,  and  Asada  [1984]), 
yielding  better  approximations  of  the  original  surface. 

4.  Set  Operations  Between  Solids 

Boundary  representation  modellers  (Braid  [1979],  Mantyla  and  Sulonen  [1982])  repre¬ 
sent  surfaces  by  faces,  edges,  and  vertices,  and  the  neighborood  relations  between  them. 
Set  operations  between  solids  are  then  computed  by  first  finding  the  boundary  intersec¬ 
tion  (which  is  a  set  of  polygons)  of  the  two  objects,  and  then  classifying  the  faces,  edges, 
and  vertices  outside  or  inside  according  to  the  operation  considered. 

The  first  step,  finding  the  intersection  polygons,  is  the  most  expensive  one.  A  naive 
implementation  leads  to  a  complexity  of  O  (m.n),  where  the  two  objects  compared  have 
respectively  m  and  n  faces.  To  speed  up  these  algorithms,  Mantyla  and  Tamminen  [1983] 
use  the  localization/proceaaing  decomposition  scheme,  by  first  localizing  the  intersecting 
faces  of  the  two  objects,  and  then  computing  the  polygons  themselves.  They  organise  the 
3D  space  for  each  polyhedron  into  a  hierarchical  data  structure,  called  Box-EXCELL. 
This  structure  is  not  intrinsic  to  the  surface.  It  consists  of  a  hierarchy  of  non  overlapping 
rectangular  cells,  each  of  them  pointing  to  all  the  faces  whose  enclosing  boxes  intersect 
it.  A  cell  that  intersects  too  many  boxes  is  subdivided,  so  that  the  number  of  boxes  by 
terminal  cells  remains  small.  For  each  edge  of  one  polyhedron,  they  find  all  the  faces 
of  the  other  polyhedron  that  it  intersects  by  visiting  the  associated  Box-EXCELL  struc¬ 
ture.  Mantyla  and  Tamminen ’s  experiments  show  that  the  localization  step  complexity 
is  reduced  to  O  (m  +  n)  (it  cannot  be  less  in  their  case  as  each  edge  of  both  polyhedra  is 
visited). 

We  use  the  same  decomposition  approach,  but  the  localization  of  the  intersection  is 
different.  Our  algorithm  is  a  direct  generalisation  of  the  Ballard’s  algorithm  (Ballard 
[1981])  for  finding  the  intersection  of  two  curves.  We  represent  the  two  polyhedra  by 
Prism  trees  and  mark  all  the  nodes  that  correspond  to  intersecting  faces  by  visiting  in 
parallel  the  two  trees  (this  is  similar  to  the  algorithms  for  Octrees  intersections  (Meagher 
[1980]),  although  the  result  obtained  by  our  method  will  not  be  approximate,  as  for 
Octrees,  but  exact).  At  each  level,  if  two  prisms  don’t  intersect,  then  the  associated 
surfaces  don’t  intersect  either,  so  the  nodes  and  all  their  descendants  are  marked  non¬ 
intersecting.  Otherwise,  if  the  two  nodes  are  leaves  and  the  associated  STb  intersect,  we 
mark  them  as  intersecting,  as  well  as  their  ancestors,  otherwise  we  also  mark  them  non¬ 
intersecting.  In  the  remaining  case  of  two  non-leaf  intersecting  prisms,  the  bigger  prism 
is  subdivided  (this  heuristic,  originally  proposed  by  Ballard  [1981]  in  the  2D  case,  has  for 
purpose  to  always  compare  prisms  of  equivalent  sites),  and  the  recursion  proceeds. 

The  processing  step  then  begins:  we  classify  the  non  intersecting  inside  or  outside 
nodes  of  one  tree  with  respect  to  the  other  one  by  classifying  one  node,  and  propagating 
its  status  to  its  connected  component  of  non-intersecting  nodes  (as  two  neighboring  non¬ 
intersecting  nodes  are  necessarily  either  both  t'nst'de  or  both  outside),  and  eliminate 
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the  inside  or  the  outside  nodes  according  to  the  operation  performed.  The  connected 
component  exploration  is  done  by  using  a  neighbor  finding  algorithm  (Faugeras  and  Ponce 
[1983,85]),  analogous  to  the  one  developed  by  Samet  [1982]  for  Quadtrees,  and  that  finds 
any  neighbor  of  a  node  at  the  same  level  of  the  tree  in  constant  average  time.  We  have 
not  yet  included  in  the  program  the  construction  of  the  intersection  polygons,  but  we 
intend  to  do  it  in  a  near  future,  by  using  a  solid  modeller  for  polyhedra  deve  loped  by 
Lanusse  [1985].  As  already  noted,  the  addition  of  this  module  will  not  change  the  overall 
complexity  of  the  method,  as  the  number  of  faces  by  leaf  is  always  bounded. 

Figure  9  shows  the  application  of  this  algorithm  to  set  operations  between  a  sphere  and 
a  double  cone,  whose  Prism  TYees  representations  contain  approximative^  1090  leaves 
each.  Tbe  computing  time  tor  the  localisation  step  is  only  25s  on  a  VAX  class  mini, 
computer.  Of  the  million  of  possible  intersections,  less  than  one  thousand  are  actually 
tested.  The  processing  step  takes  about  Is. 

We  may  make  some  remarks  about  the  complexity  of  the  localisation  procedure.  At 
a  given  level  of  the  tree,  only  the  intersecting  nodes  are  again  subdivided.  This  would 
lead,  if  only  one  of  the  sons  of  each  intersecting  node  was  also  an  intersecting  one,  to  a 
complexify  of  O  where  TV,-  is  the  number  of  intersecting  faces,  and  K  the  depth 

of  the  tree.  Such  a  best  case  is  evidently  not  realised  in  practical  experiments,  but  it 
indicates  that  the  complexity  depends  on  the  shape  of  the  intersection  polygon,  and  on 
the  logarithm  of  the  total  number  of  faces  (corresponding  to  the  depth  K  of  'he  tree), 
so  the  method  seems  better  than  the  Mantyla  and  Tamminen  one,  although  a  more 
careful  study  is  necessary.  Notice  anyway  that  an  advantage  of  this  method  is  that, 
the  representation  being  intrinsic,  it  does  not  have  to  be  recomputed  when  the  object  is 
rotated  (the  same  advantage  holds  with  respect  to  Octrees  (Iftikhar  [1981]).  A  procedural 
version  of  the  localization  algorithm  is  given  in  the  appendix. 

5.  Ray  Casting 


Hidden  surface  elimination  algorithms  (Sutherland,  Sproull,  and  Schumacker  [1974]), 
and  in  particular  Ray  Casting  methods  (Roth  [1982]),  although  they  are  among  the  more 
general  and  allow  a  lot  of  “special  effects” ,  are  usually  computationally  expensive:  if  Ns  is 
the  number  of  surface  elements,  Nt  ray-surface  intersections  must  be  computed  for  each 
pixel.  This  makes  the  “naive”  version  of  this  approach  unusable  for  moderate!)’  complex 
polyhedral  objects  each  often  containing  more  than  1000  facets  for  images  containing 
at  least  a  quarter  million  pixels.  To  increase  the  efficiency  of  Ray  Casting  algorithms, 
localization  methods  have  been  used. 

Clark  [1976],  and  after  him  Rubin  and  Whitted  [1980],  enclose  the  objects  composing 
the  scene  in  a  hierarchy  of  simple  boxes  (parallelepipeds).  At  each  level,  if  the  box 
does  not  intersect  the  ray,  the  corresponding  objects  are  rejected.  Otherwise,  the  box  is 
subdivided,  etc..  Weghorst,  Hooper,  and  Greenberg  [1985]  have  successfully  applied  more 
sophisticated  versions  of  similar  methods  to  very  complex  scenes.  However,  these  methods 
are  object  oriented:  they  use  hierarchies  of  objects,  bounded  by  simple  volumes,  but  they 
are  not  really  designed  for  dealing  with  polyhedra  that  contain  themselves  thousands  of 
faces. 


Figure  10.  The  fractal  subdivision  proccse:  for  eoch  edge  of  the  original  triangle  T,  a  subdivision  point  of 
height  Mi  is  computed,  and  T  is  replaced  by  the  four  triangles  2ms,  Tu,  Tit,  and  2||.  The  “cheesecake’’ 
is  obtained  by  translating  the  original  triangle  along  the  Z  aaris. _ 

Recently,  Kajiya  [1983]  has  proposed  a  face  oriented  method  for  rendering  fractal 
surfaces  (Fournier,  Russel  and  Carpenter  [1982])  using  Ray  Casting.  His  surface  model 
is  similar  to  ours,  except  that  it  is  non  deterministic.  The  surface  is  a  polygonal  height 
field.  It  is  recursively  subdivided,  each  triangle  being  replaced  by  four  subtriangles  whose 
vertices  are  generated  by  a  stochastic  process  (the  subdivision  in  4  triangles  eliminates  the 
need  for  an  adjustment  step,  but  on  the  other  hand,  the  subdivision  cannot  be  stopped 
for  some  triangles  as  the  polyhedral  structure  would  not  be  preserved).  Kajiya  encloses 
each  triangle  in  a  so  called  “cheesecake”  extent  (Figure  10),  which  is  a  prism  obtained  by 
translating  the  triangle  along  the  vertical  to  enclose  the  associated  surface.  To  intersect 
the  ray  with  a  surface,  Kajiya  visits  the  associated  tree:  if  the  ray  intersects  a  box,  then  it 
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is  again  tested  against  tbs  descendants.  Otherwise,  the  ray  is  guaranteed  not  to  intersect 
any  of  the  descendants.  The  tree  is  visited  until  the  first  intersecting  leaf  triangle  is 

The  method  is  readily  extendible  to  Prim  7 bees.  There  are,  however,  some  differ¬ 
ences:  as  the  fractal  surface  is  defined  as  a  height  field,  the  cheesecakes  associated  to  the 
sens  of  a  given  node  never  intersect  each  other.  This  implies  that  they  can  always  be 
sorted  from  the  nearest  to  the  farthest  with  respect  to  a  given  ray.  In  particular,  if  only 
one  fractal  surface  is  displayed,  the  corresponding  tree  can  be  visited  in  a  strict  depth 
first  manner.  If  the  closest  node  is  always  visited  first,  then  it  is  guaranteed  that  the  first 
intersection  found  will  be  the  actually  closest  one. 

Prim  Trees  do  not  have  this  nice  property,  and  in  particular  non-regular  sons  of  a 
same  node  may  intersect.  We  have  to  maintain  an  active  list  of  all  the  nodes  whose 
associated  surface  may  intersect  a  given  ray.  This  disadvantage  is  a  price  to  pay  for 
generality.  Notice  however  that  the  same  problem  would  arise  for  fractal  surfaces  defined 
on  non-planar  objects,  or  simply  when  several  of  them  are  to  be  displayed  simultaneously. 
Moreover,  the  clear  intersection  lemma  of  Section  3  is  going  to  help  us  keeping  the  active 
list  short. 

Assume  that  the  ray  is  parameterised  by  A,  and  let  Xmn  ( Pt )  and  A**,!  [Pt)  be  the 
values  of  this  parameter  at  the  extremities  of  the  intersection  of  the  ray  with  the  node  Pt. 
Following  Kajiya,  we  say  that  a  node  Pt\  shadows  a  node  Pt%  if  Amat  (Pti)  <  Aran  (Pt *) 
(Figure  11).  Kajiya  conjectures  that  this  notion  could  be  used  to  prune  the  active  list, 
but  remarks  that,  unfortunately,  a  node  can  be  shadowed  at  a  given  level  and  though 
have  visible  descendants,  as  the  intersection  of  the  ray  with  the  box  does  not  imply  the 
intersection  of  the  ray  and  the  surface.  We  solve  this  problem  by  proving  the  following 
lemma. 

Shadowing  Lemma:  if  the  intersection  of  Rag  and  Pt\  is  clear ,  and  if  Pti  is  shadowed 
by  Pti,  then  Ptj  and  his  descendants  are  not  visible  from  the  considered  pixel,  and  Ptj 
can  be  removed  from  the  active  Hat. 

Proof:  as  the  intersection  is  clear,  Ray  intersects  the  surface  associated  to  Pti  at  a 
point  that  verifies  A,^  (P*i)  <  A  <  Amax  (Pti),  and  so  hides  the  surface  associated  to 
Pti,  whose  each  point  verifies  Amu  (Pti)  <  A***  (Ptj)  <  A  a 

Using  this  property,  we  can  find  for  each  pixel  the  point  of  the  scene  seen  by  the  pixel  by 
a  breadth  first  visit  of  the  trees  of  the  objects  composing  the  scene  (so  as  to  eliminate  as 
early  as  possible  the  shadowed  nodes).  At  each  step  of  the  recursion  a  A™,* -sorted  list  of 
the  nodes  that  intersect  the  ray  at  this  level  is  maintained.  The  process  terminates  when 
the  list  is  either  empty  or  only  composed  of  leaves.  The  first  element  of  the  list  is  then 
used  to  compute  the  display  parameters,  by  finding  among  the  (few)  faces  associated  to 
this  leaf  the  closest  one  that  intersects  the  ray. 

Again,  we  discuss  the  complexity  of  this  algorithm.  First  notice  the  interest  of  the 
clear  intersection  notion.  In  the  ideal  case,  there  are  only  clear  and  null  intersections,  so 
the  algorithm  visits  only  one  branch  per  tree  and  per  pixel,  and  due  to  the  shadowing 
between  trees,  one  can  expect  to  visit  entirely  a  single  branch  of  a  single  tree.  This  is  an 
important  advantage  over  K^jiya’s  algorithm,  where  the  execution  time  is  proportional 


to  the  number  of  objects  in  the  scene  (although  it  is  unlikely  that  several  fractal  surfaces 
are  going  to  be  displayed  in  the  same  image). 

In  the  case  where  few  clear  intersections  occur,  or  the  different  surfaces  are  close 
enough  so  that  only  few  shadowings  occur,  the  complexity  degenerates,  and  becomes 
again  proportional  to  the  number  of  objects.  Even  in  this  case  though,  it  is  likely  that  a 
ray  will  in  general  intersect  only  one  of  the  sons  of  a  node,  so  only  a  few  branches  will  be 
explored. 

Figure  12  shows  the  application  of  this  algorithm  to  the  union  of  the  double  cone 
and  the  sphere,  and  to  a  Renault  automobile  part.  The  resolutions  of  the  images  are 
respectively  512  x  512,  with  a  CPU  time  of  two  hours  for  2000  leaves,  and  256  x  256,  with 
a  CPU  time  of  28  minutes  for  1000  leaves.  Once  again,  in  our  case,  the  processing  step  is 
not  complete,  as  we  display  the  approximation  and  not  the  surface  itself.  Nevertheless,  the 
remark  made  at  the  end  of  the  previous  section  still  holds,  and  a  complete  implementation 
will  not  change  the  overall  complexity  of  the  algorithm.  A  procedural  version  of  the 
algorithm  can  be  found  in  the  appendix. 

0.  Conclusion 

We  have  presented  a  new  method  for  localizing  the  search  for  geometric  intersections 
in  the  polyhedral  case.  It  has  proved  efficient  for  the  two  classical  problems  studied.  We 
are  continuing  our  work  on  the  complete  implementation  of  the  processing  step.  Also,  a 
more  rigorous  study  of  the  algorithms  complexity  is  needed. 
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Figure  12.  Examples  of  applications  of  Ray  Casting  to  the  union  of  the  double  cone  and  the  sphere,  and 
to  an  automobile  part 


We  believe  that  Prism  Trees  have  important  potential  applications,  in  particular  for 
complex,  real  3D  objects,  digitized  as  polyhedra  with  thousands  of  facesf.  In  this  optic, 
we  give  two  examples. 

Today’s  most  sophisticated  Graphics  systems  deal  with  (relatively)  simple  objects.  Us¬ 
ing  the  Ray  Casting  algorithm  in  conjonction  with  a  laser  rangefinding  system  (Faugeras 
and  Pauchon  J983  )  will  allow  us  to  include  “real  life”  objects,  from  industrial  parts  to 
statues,  in  such  an  high  level  Graphics  environment  and  to  obtain  from  them  the  realistic 
images  that  only  Ray  Casting  programs  authorize. 

Similarly,  the  Set  Operations  algorithm  could  be  used  in  CT  tomography  (Artzy, 
Frieder  and  Hermann  j  19781)  to  obtain  slices  of  an  organ  in  an  arbitrary  direction,  or 
more  generally  to  cut  it  along  an  arbitrary  surface. 

7.  Appendix 

We  now  present  in  more  detail  the  procedures  used  to  find  the  intersection  of  two 
surfaces,  and  to  display  the  object  represented  by  the  Prism  Tree  using  a  Ray  Casting 
method.  We  give  them  in  a  pseudo  PASCAL  form,  the  derlarations  of  types  and  variables 
being  omitted. 

7.1.  The  Surface-Surface  Intersection  Procedure 

This  procedure  is  straightforward,  the  two  trees  are  visited  in  parallel  in  a  depth 
first  search  manner,  until  all  the  intersection  leaves  tire  marked  (using  the  flag  Mark). 


f  This  is  the  typical  output  of  scanning  processes  that  measure  successive  slices  of  the  surface  of  a  3D 
object,  as  laser  digitizers  or  CT  scanning. 
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Then  ail  their  ancestors  are  also  marked.  We  suppose  that  we  dispose  of  the  following 
functions; 

Vol  [Pt):  this  real  function  returns  the  volume  of  the  boot  associated  to  Pt. 

Leaf  (Pt):  this  boolean  function  returns  true  iff  the  node  Pt  is  a  leaf. 

Marked  (Pt):  this  boolean  function  returns  true  iff  Pt  is  not  ntf  and  is  marked  (i.e.  the 
flag  Mark  is  set  to  1). 

Priemint  (PthPtj)  is  a  3*  valued  function  that  returns  clear,  possible,  or  null  according 
to  the  type  of  intersection  of  Pt\  and  Ptj.  This  function  is  the  key  of  the  algorithms  that 
manipulate  the  Prism  TVees.  In  the  special  case  of  the  leaves  whose  associated  prisms 
{  intersect  (and  only  in  this  case),  Prismint  tests  directly  the  intersection  of  the  surface 

patches  themselves,  and  returns  clear  if  they  actually  intersect. 

Procedure  Find  -  Surfaces  -  Intersection 

V - i— 

Dcgm 

if  Prismint  (Ptt,Pt,)  ?  null  then 

h»gtn 

Lei  Leaf  (Pti);  Le,  «-  Leaf  ( Pt,); 
it  Lei  and  Let 

then  begin  Pt,  f  .Mark  *-  1;  Pt,  |  Mark  *-  1  and 
else  (*  The  recursion  proceed*  *) 

if  Le,  or  ( Vol  (Pt,)  >  Vol  (Pt,)) 
then  begin 

for  t «-  1  to  3  do 

Find  -  Surfaces  -  Intersection  (Pt,  |  .Son  [t]  ,Pt,); 
if  Marked  (Pt  , 1  .Son  [1]) 
or  Marked  ( Ptj  f  .Son  [3]) 
or  Marked  (Pt  j  f  .Son  [3]) 
then  Pti  t  Mark  «-*  1 
end 
dae  begin 

for  < «-  1  to  3  do 

Find  -  Surfaces  -  Intersection  (Ptj,Pt,  t  .Son (»']); 
if  Marked  (Pt,  f  .Son  [l]) 
or  Marked  (Pt,  f  .Son  [2]) 
or  Marked  (Pt,  t  • Son  [3]) 
then  Pt,  t  Mark  1 
end 
end 

end; 


7.2.  The  Ray  Casting  Procedure 

The  elements  of  the  Amtn-sorted  list  of  nodes,  Active  -  List,  ere  5*fields  records:  the 
first  field,  Node,  is  the  corresponding  node  of  the  Prism  Tree,  the  second  one,  Next,  is 
a  pointer  to  the  next  element  of  the  list,  the  third  one,  Old  is  used  for  the  breadth  first 
visit  of  the  tree.  Its  value  is  1  if  the  node  has  been  inserted  at  the  previous  iteration. 
The  last  two  fields  caracterise  the  intersection  of  the  ray  with  the  prism.  X^ode  >s  the 


V  *  » 

i**  *  . 


■  •  '  «  •  •  •  • 


Prism  Trees 


minimum  nhic  of  the  parameter  A  of  the  ray  each  that  an  intersection  occurs,  and  Pfiodei 
set  only  if  the  node  is  a  leaf  and  has  a  clear  intersection  with  the  ray,  is  the  firstpoint  of 
intersection  of  ST  and  the  ray.  It  is  used  for  computing  the  display  parameters. 

We  suppose  to  dispose  of  the  following  procedure  and  function: 

Insert  (Q,JV,Afm«)  in*«rts  the  element  Np  in  the  list  Q  by  comparing  the  A***  value 
to  N  to  the  one  of  the  currently  visited  node. 

Rayint  (Ray,  Pt ,  A***,  Xmas,  Pmin) »  a  modified  version  of  Prismint .  It  is  also  a  3-vafoed 
function,  the  clear  intersection  notion  being  extended  to  the  intersection  of  a  straight  line 
anri  ^  regular  prism.  Its  additional  features  are  that  it  returns  the  extremal  values  A^o 
and  of  the  intersection  of  Ray  and  Pt,  and  in  the  case  where  Pt  is  a  leaf  and  the 
intersection  is  clear ,  it  also  returns  the  closest  intersection  point  P***. 

Function  Rayeast  (Ray,  Pt,  P,*«n)  boolean; 
begin  (*  initialize  parameter*  *) 
with  Active  —  List  T  .Next  f  do 

begin  Node  «—  Pt;  Old  «—  1;  A jvo*  end; 

Am«x,o  +oo;  A  4  +oo;  Ptttn 
repeat  (*  main  loop  *) 

All  -  Leave*  true;  Empty  «-  true; 

P  «—  Active  -  List ;  Q  *-  P  t  .Next; 

while  Q  ft  nil  do  (*  auxiliary  loop:  visit  the  list  *) 

if  Q  T  AjVWe  >  Am,,0 

then  Q  «—  nil  (*  the  rest  of  the  list  is  shadowed  *) 
rise  if  Leaf  (Q  t  .Node)  (*  non  shadowed  leaf?  if  yes,  it  shadows  *) 

then  begin  Empty  <—  false;  Q  <—  nil  end  (*  everything  behind  *) 
else  if  Q  T  • Old  =  0 

then  (*  non-leaf  new  node  that  gets  old,  go  on  *) 
begin  All  -  Leaves  <—  false;  Q  T  . Old  —  1; 

P  Qi  Q  Q  T  .Next  end 
else  begin  (*  non  leaf  old  node,  divide  *) 
for  t  1  to  3  do 
begin 

Inter  «-  Rayint  (Ray,  Q  T . Node  ]  .Son  [«]  ,Am<n,  Amoxi  P mi* 
if  ( Inter  #  null)  and  (Am,„  <  Amax.o)  then 
begin  (*  insert  the  son  *) 

New(Np);  Insert  (Q,Np,\min)\ 
with  Np  T  do 

begin  Old  *—  0;  A fiodt 4  Am,„;  Pno* u  P mini 

Node  Q  t  .Node  t  . Son  [i]  end; 
if  ( Inter  =  clear)  and  (Am„  <  Am.*,o) 
then  Am,,,!) 4  Amal ; 

end  end; 

Q  «-  Q  T  .Next;  P  f  .Next  *-  Q 

end; 

until  All  -  Leaves  at  Empty; 
if  Empty  then  Rayeast «—  false 

else  begin  PMtn  Active  -  List  |  .Next  |  J*NoSt\  Rayeast  «—  true  end 
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